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Abstract 


Experimental and theoretical results are presented for nonlinear dispersive 
waveguide modes in an acoustic duct. We report observations of a localized 


envelope which Splits into two disturbances moving with two different velocities 


of propagation. As a consequence, we predict that if a signal is amplitude 


modulated at the source, spatial beating between the two disturbances will 
occur, and at periodic positions in space the signal will become frequency 
modulated. These results can have applications in an all-optical AM-FM 


conversion and in high data rate fiber optic communications. 
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CHAPTER 1. INTRODUCTION 


Nonlinearity can give rise to behavior that is fundamentally different from 
linear behavior. Solitons are the most dramatic example where an exponentially 
localized wave of constant shape is the result of a stable balance between 
nonlinearity, which causes the wave to shock, and dispersion, which causes the 
wave to spread. Among the most well known examples of solitons are the 
envelope solitons. These solitons are the result of an instability whereby an 
amplitude modulation superimposed on a wavetrain initially grows and eventually 
self-localizes. As we will see below, this instability is realized if the product of the 
nonlinear coefficient and the curvature of the dispersion curve (frequency vs. 
wavenumber) is negative. 

When the product of the nonlinear coefficient and the curvature of the 
dispersion curve is positive, modulations of a nonlinear wave are stable. In this 
case, the most striking feature resulting from the nonlinear action is that a 
modulation of finite extent will split into two disturbances moving with two 
different velocities of propagation (Whitham, 1974). This behavior should be 
contrasted with linear theory where such a disturbance will undergo distortion 
due to dispersion, but will not split. To our knowledge, no experimental evidence 
of this effect has been pursued to the present, with the exception of this work. 

We can gain a physical understanding of amplitude modulations of nonlinear 


dispersive waves by considering independently the effects of dispersion and 





nonlinearity on a wave. Although the ideas below are very general, we will 
develop the conceptual background for acoustics in a waveguide, as this is the 
system that we use for our experiment. 

A harmonic plane wave propagating in free space or in a waveguide below 
cutoff satisfies the dispersion relation m = ck , where @ is the frequency of the 
wave, k is the wavenumber of the wave, and c is the speed of sound in the 
medium. That is the frequency is proportional to the wavenumber and the phase 
velocity equals the slope dw/dk of the line in the w-k plane. In this case, the 
phase and group velocities are independent of frequency. Therefore, a group or 
superposition of harmonic waves with different frequencies all travel together 
with the same speed. 

However, for acoustic propagation in a waveguide above cutoff the dispersion 
relation has curvature (i.e., d?w/dk? =a” (k) #0). Hence the phase velocity is 
frequency dependent and a group of waves of different frequency - and therefore 
speeds - will change shape as it propagates. Since the motion of the wave 
group can no longer be described by a single velocity, a group velocity v,=dw/dk 
is defined. The group velocity gives the velocity of propagation of a group of 
waves extending over a frequency range small compared to the center frequency 
of the group. 

To clarify the concepts of the phase and group velocities, consider an 
observer in a frame moving with the phase velocity. In this frame an 
unmodulated wave is described by a stationary pattern of crests and troughs 
fixed in space and the different points on the wave correspond simply to points of 


different phase (Fig. 1.1a). To an observer in this frame, a modulated wave 








would also appear as a pattern of crests and troughs fixed in space (Fig. 1.1b). 
However if w’’(k) # 0, energy flows from regions of high to low amplitude with 
the group velocity and the amplitude of the crests and troughs would vary in 


time. Thus an observer fixed in a frame moving with the phase velocity would 


observe a propagating modulation with a relative speed |v,-v,]|. 








a) 


ue 


b ) 


\ energy density 








a energy density 


Fig. 1.1 a) Unmodulated and b) modulated wave in a frame moving 
with the phase velocity v,. 


For large amplitude acoustic waves nonlinear effects become important. For 
plane waves in a nondispersr. 2 medium, distortion of the sound wave results 
from amplitude dependent effects. Since the wave speed is proportional to the 
square root of the temperature, it follows that compressions, where the 
temperatu.e is higher, propagate faster than rarefactions, where the temperature 
is lower. Also, at the compression the particle velocity is in the direction of 
propagation and at a rarefaction the particle velocity is in the opposite direction. 
Shock formation of the waveform results in the sawtooth shape shown in Fig. 
1.2a as the crests of the wave overtake the troughs. 

For a dispersive medium shocks are not produced because frequency and 
wavenumber are not proportional and thus new propagating modes are not 
generated by nonlinear interactions. However, due to self interaction effects, the 
frequency of the wave is amplitude dependent and in the weakly nonlinear 
regime it is of the form @ = w,(k) + «(k)p?. In this expression, o,(k) is the linear 
dispersion relation, w,(k) is the nonlinear coefficient, and p is acoustic pressure. 
For tlie case of fixed frequency and w, > 0, the effect of nonlinearity is to 
decrease the wavenumber k as the amplitude of the acoustic field increases. 
Equivalently, for fixed frequency, as the amplitude of the wave is increased the 
phase velocity increases. 

Because the changes in the wavenumber vary linearly with p’, the relative 
phase between two points fixed in space should also vary in this manner. In 
Chapter 5 we present experimental results that verify this effect for the (1,0) 
cylindrical waveguide mode. This result is the first experimental observation of 


the dependence of the relative phase on the wave amplitude. 


a) 
before 
after 
b) 
© = @,(k) + @2(k)p? @>0 


Fig. 1.2 a) Shock formation resulting from nonlinearity in a 
nondispersive medium as viewed in a frame moving with 
the speed of sound. b) For w,>0 and fixed frequency, the 
effect of nonlinearity in dispersive acoustics is to increase 
the wavelength. 





To understand the combined effects of dispersion and nonlinearity, consider 
an initial state of a modulated wave observed in a frame moving with the group 
velocity (Fig. 1.3a). Because of dispersion, the group and phase velocities are 
different and an observer in this frame would see the crests of the wave 
propagate. If the phase velocity is bigger than the group veiocity (as is the case 
for a mode above cutoff in a waveguide) the crests of the wave will propagate in 
the positive x direction. 

Without loss of generality we can assume that the nonlinear coefficient w,>0. 
Thus, the phase velocity in the regions of higher amplitude is greater than in 
regions of lower amplitude. As a result, an observer in the frame moving with the 
linear group velocity would observe bunching of the crests at the leading edge of 
the modulation and anti-bunching at the trailing edge (Fig. 1.3b). For negative 
dispersion, w"(k)<0, the long waves travel faster than the short waves and an 
instability sets in. Dispersive effects cause the energy to approach the crests of 
the envelope, and the amplitude at the crest increases. This in turn accelerates 
the instability. Deep gravity waves on the surface of a liquid is an example of a 
wave system with negative dispersion. For positive dispersion, w’(k)>0, short 
waves travel faster than long waves. In this case dispersive effects cause the 
energy to approach the troughs of the envelope and the modulation propagates. 
In other words, to an observer moving with the linear group velocity, the 
modulation is no longer stationary. Acoustic propagation in a waveguide above 


the cutoff frequency is an example of a wave system with positive dispersion. 


In general, if the product ,” is negative a modulation will initially grow and 


the instability will eventually lead to the formation of envelope solitons. {f the 





product ww" is positive, the modulations are stable. In Chapter 4 we calculate 
the nonlinear coefficient for the (1,0) mode in a cylindrical waveguide and show 
that it changes sign from positive to negative as the frequency is increased. For 
this mode, finite extent modulations of a uniform wavetrain are stable below a 
frequency approximately twice the cutoff frequency and are unstable beyond this 


frequency . 


— - 


Fig. 1.3 a) Initial state of a modulated wave in the frame moving 
with the group velocity v,. b) For w,>0, in a frame moving 
with v, nonlinearity results in bunching of the crests at the 
leading edge of the modulation and anti-bunching of crests 
at the trailing edge. 





The stability of the modulation has two important consequences. One is 
AM-FM conversion. The other is pulse splitting. Both effects can be physically 
understood by considering again a modulated wave in a frame moving with the 
linear group velocity. !n this frame, an observer at a fixed location would see 
alternating bunching and anti-bunching of the wave crests (Fig. 1.4). Because 
dispersive effects cause the energy to approach the troughs of the envelope, for 
this observer it would appear that at some time a frequency modulation is 
imposed upon the amplitude modulated signal. The combined dispersion and 
nonlinearity cause the energy flow from the short wavelength to the long 
wavelength of the FM modulation and an amplitude modulation 180° out of 
phase with respect to the original signal results at a later time. In the frame 
moving with the group velocity, the process repeats periodically and an observer 
in this frame sees that the modulations experience beats. In the laboratory 
frame if a source is generating an amplitude modulated signal, some distance 
away it will become frequency modulated. This spatial beating of the modulation 
can be only be understood if for each modulation frequency there correspond 
two group velocities of propagation. 

The splitting of a single high intensity modulated wave has been verified 
experimentally and is reported in Chapter 5. A single pulse modulation on a CW 
is observed to split into two pulses. The time difference between them increases 


both with distance and amplitude. 
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Time 





Fig. 1.4 An observer in a frame moving with the linear group 
velocity observes AM to FM conversion in time. 
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The organization of this work is as follows. In Chapter 2 we derive the 
coupled wavenumber and energy conservation equations. Chapter 3 deals with 
the quantitative derivation of group velocity splitting. As a consequence of group 
velocity splitting we predict the existence of AM-FM conversion in modulationally 
stable nonlinear dispersive systems. In Chapter 4 we apply these results to the 
specific case of the (1,0) mode in a waveguide of circular cross section. A 
calculation of the nonlinear coefficient w, for this mode is presented and the 
experimental consequences of these results are analyzed. In Chapter 5 we 
present the experimental results for a nonlinear (1,0) mode in an acoustic duct 
for the frequency regime that leads to modulational stability. Finally, Chapter 6 
summarizes the results of our theoretical and experimental investigations of this 


phenomena. 
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CHAPTER 2. ONE DIMENSIONAL MODULATION THEORY 


in this chapter we review the theory of one dimensional modulation of linear 
waves in a dispersive medium. We then expand our discussion to include the 
effects of nonlinearity and the combined effects of both nonlinearity and 


dispersion. Our treatment parallels the work done by Whitham (1974). 


1. LINEAR MODULATION THEORY IN ONE DIMENSION 
a) Kinematic Considerations 
If the amplitude and propagation wavenumber of a traveling wave vary 
only slightly over distances of the order of one wavelength, the wave can be 
described locally by assuming a slowly varying amplitude and phase. A wave 


field ‘¥Y can then be described by the expression 


P=a(x,the (2.1.1) 


where a(x,t) is the local value of the amplitude. Over small regions of space and 
short intervals of time the phase (x,t) can be expanded in a Taylor series and 


up to first order, we obtain 


ox at 


where 6, is the phase evaluated at the origin at time t=0. 





Over small regions of space and during short time intervals, we define the 


local wave number and frequency at each point as 
o,5) 6,8) 
k(x, t) = — ; v(x, t)=-—. 2.1.3 
(x,t) x (x,t) a ( ) 
Thus, points of constant phase move with velocity 
ee (2.1.4) 


An observer following a particular crest moves with the local phase velocity but in 
general sees the local wave number and frequency changing; that is neighboring 
crests get farther away. 


By eiiminating the phase 6 from each expression in Eq. (2.1.3), we obtain 
a a (2.1.5) 


Equation (2.1.5) has the form of a conservation law, where the 
wavenumber k in the first term corresponds to the density of waves or number of 
wave crests per unit length and the frequency in the second term corresponds 
to the wave flux or number of wave crests crossing the position x per unit time. 

Equations (2.1.1) - (2.1.3) can be extended to include the case of two and 
three dimensions. However, the one dimensional equations are sufficient to 
describe waveguide propagation and we will use them exclusively throughout our 


discussion. 
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b) Energy Propagation 
For a general dispersion relation w = w(k) the phase velocity is in general 
frequency dependent. Thus, if a disturbance consists of a group of individual 
harmonic waves of different frequencies and hence speeds, the resulting wave 
shape will change with position and the motion of the wave group cannot be 
described in terms of a single wave velocity. 
The general form for a propagating wave in a linear dispersive medium 


with dispersion relation w(k) is 
(x,t) = [F(k)eOM'dk (2.1.6) 


For a modulated wavetrain with most of the energy in wavenumbers close to 
some value k_. F(k) is concentrated near k=k, and (2.1.6) may be approximated 


by 
(x,t) = JF (k)exp|[(ikx -i{o(k,) + (kk, )o'(k,)}t]dk . (2.1.7) 
This in turn may be written as Y = a(x -0'(k, )t)exp{i(k,x—o(k,)t} where 
a(x—o'(k,)t) = [F(k, +)explix(x-a'(k,)t}dx (2.1.8) 


and where we have substituted k = k, + x. The function a(x—'(k,)t) describes 
the modulations of the wave field about k,. It varies with x and t in such a way 
that it remains constant when viewed from a coordinate system moving with the 
group velocity 


_ da 


Lae 2.1.9 
Ca ak ( ) 
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evaluated at the center freyuency of the group, i. e. at @ =w(k,). Thus, we say 


that the modulation is constant along the group line x=o'(k,) t. An observer 
moving at the group velocity w'(k,) always sees waves with wavenumber k, and 
frequency w(k,) but crests keep passing him. 

The energy density of a wave is proportional to the square of the wave 
amplitude. Hence if a(x—w’'(k, )t) is constant in a coordinate system moving 
with the group velocity, the energy between any two group lines remains 


constant. In other words, 


ce oO 

—+—(c,e)=0 , 2.1.10 

is a g@) ( ) 
where e is the energy density. In order to see that (2.1.10) is valid, consider the 


quantity 


E(t) = in a’dx = 2n [CF QOF(k) ‘dk, (2.1.11) 


where k, and k, are defined by 
xX, =@'(k,)t ’ Xo ='(kK,)t . (2.1.12) 


If k, and k, are held fixed as t varies, E(t) remains constant. According to 
(2.1.12), the points x, and x, are moving with the corresponding group velocities. 
Thus the total amount of a” between any pair of group lines remains the same. 
Thus a? propagates with the group velocity. On the other hand, e = f(k) a”, thus 
(2.1.10) holds provided 








da 8 ay | ga yn2fOkK , . Ok 
10 + Soya yp (k)a {B+e, Xf =o. (2.1.13) 
The first bracket vanishes because E(t) is constant between any pair of group 
lines. The second bracket vanishes in view of (2.1.5), the conservation of the 


number of waves crests! 


2. NONLINEAR WAVES 

The modulations of a linear wavetrain propagating in one dimension in a 
uniform medium can be described by Eqs. (2.1.5) & (2.1.10). More generally, 
nonlinearities modify the wave propagation. Nonlinearity results in the 
appearance of an amplitude dependence ‘= the frequency. For small 


amplitudes, the frequency dependence on amplitude can be written as 
w(k,a?) =@9(k)+@,(k)a*+... , (2.2.1) 


where a is the amplitude, , (k) is the linear dispersion relation, and w,(k) is the 
nonlinear coefficient. Since the energy e of the wavetrain is proportional to the 


amplitude squared, Eq. (2.2.1) can also be written 
wo(k,e) = 09 (k)+@>2(k)e+... . (2.2.2) 


There are three conditions possible for the nonlinear coefficient w, . If 


> (kK) =0, the dispersion relation reduces to the linear case. !f m2(k)< 0 both 
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the frequency and phase velocity decrease with increasing amplitude. Finally, if 
(kK) > 0, the frequency and phase velocity increase with increasing amplitude. 
An important effect of the nonlinearity is to couple Eqs. (2.1.5) and (2.1.10). 
Thus substituting (2.2.2) into (2.1.5) we obtain 
ok 


ok de 
St (ook) + a2 (k)e) + @9(k) = = 0. (2.2.3) 


For slow modulations of the amplitude, the last term of Eq. (2.2.3) can be as 
significant as the term «4 (kK) Ok/Ox. However, the w,(k)e0k/dx term can be 
ignored for weakly nonlinear fields because it only gives a small correction to the 
group velocity. 

Similarly, in Eq. (2.1.10), we can retain only the linear group velocity and 


obtain the equations (Whitham, 1974) 


dk ék de 

at +@o(k) = + o2(k) = = 0 : (2.2.4) 
oe oe ck 

a POalk) = +9 (kK)e= = 0. (2.2.5) 


Equations (2.2.4) and (2.2.5) constitute the basic description for one dimensional 
nonlinear waves in a uniform dispersive medium. The crucial qualitative change 
of nonlinearity which couples both equations will be explained in the next 


chapter. 
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CHAPTER 3. GROUP VELOCITY SPLITTING 


In this chapter, we consider the general solutions to Eqs. (2.2.4) and (2.2.5). 
Based on Whitham's modulation theory we will give a quantitative foundation for 
the splitting of a modulation of finite extent into two separate disturbances. Next, 
in an extension of these results, we solve the eigenvalue problem for a 
monochromatic space-time modulation of wavenumber and amplitude. From this 
we determine that the modulational energy is equipartitioned into two 
propagating modes. Lastly, we predict that a consequence of nonlinearity on an 
amplitude modulated signal is to introduce a frequency modulation into the signal 


and that the reciprocal effect is also possible. 


1. PULSE SPLITTING 

The combined effects of dispersion and noniinearity on an initial disturbance 
or a modulation is to couple Eqs. (2.2.4) and (2.2.5) 

This system of equations is hyperbolic and in order to determine the 
characteristic velocity for the system (2.2.4) - (2.2.5), i.e. the velocity on the 
characteristic along which the wave propagates, we consider the linear 


combination® 


& For an excellent treatment on the classification of partial differentia! equations and the 
determination of the family of characteristic curves for hyperbolic equations, the reader must 
consult Whitham's Linear and nonlinear waves . 








Ok ek oe oe ce ok 
4{ S+05) E+ 0,0) S }o1,[ F roo tky S +09 (kd )= 0, (3.1.1) 


ot ot 


which takes the characteristic form 


(+E E) (SRE I-0, (3.1.2) 
at dt ox at dt &x 
provided 
I, eae 14 (k) +loog (ke, 
. (3.1.3) 
x 
bee = Ho (k) +hera(k). 


The system of linear algebraic equations (3.1.3) possesses a non-trivial solution 
when 
dx 


o% = oh (k) +(w2(k)o4 (k)e) 


% 
dt 


(3.1.4) 


and the solutions are real (i.e. propagation is possible) provided w,(k)o{ (k) 20. 
Tnus an initial disturbance or a modulation of a source will introduce 
disturbances into both families of characteristics and consequently two signals 
moving with different velocities will be detected. This behavior should be 
contrasted with the linear theory where such a disturbance will undergo distortion 


due to dispersion, but will not split. 
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2. THE EIGENVALUE PROBLEM 

In the previous section we determined that for o,(k)o ; (k) 2 0, two real 
propagating solutions to Eqs. (3.1.3) exist. In this section, we approach the 
solution of the coupled equations (2.2.4) - (2.2.5) by considering the condition on 
the values of possible frequencies. 

Consider a wavetrain which has been given a space-time modulation of 


wavenumber and amplitude: 


k(x,t) = ky +a(x,t) , (3.2.1) 
e(x,t) =e, +£(X,t) , (3.2.2) 


where k, is the unmodulated wavenumber, a(x,t) is the modulation of the 
wavenumber, e, is the constant energy of the unmodulated carrier, and (x,t) is 
the energy of the modulation. 

A linearization of Eqs. (2.2.4) and (2.2.5) about the uniform state k, and e, 


yields the modulation equations 


0a ; Ca Of 
+ 


at eta a a ’ (3.2.3) 
6:3 Oe 0a 
F100 ay tue z= 0. (3.2.4) 


if we seek solutions to Eqs. (3.2.3) - (3.2.4) of the form a =a,e"™ and 


e=e,e*"™ | we obtain 
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(Q-o)K)ag-o Keg =0, (3.2.5) 


-0f pKa +(Q-ahxe, =0 (3.2.6) 
Equations (3.2.5) - (3.2.6) have a non-trivial solution only if 


Q=0/K + ojO,e K. (3.2.7) 


The dispersion relation (3.2.7) shows that for fixed Q , there are two possible 


wavelengths of the modulated state i.e. 


Q 


q= ———_———. , 
(of, + for/er.€, ) 


(3.2.8) 


Q 


= . 
ee 
(ons ofor€o) 


(3.2.9) 


Because the nonlinearity is weak (i.e., (@9@ €) <<), (3.2.8) and (3.2.9) 


can be rewritten as 


if 
S{-seloe | (3.2.10) 


ar rs 
i 
p= 8 [Leto (3.2.11) 
@ @ 
0 0 


22 


Because there are two eigenvalues to the set of Eqs. (3.2.3) - (3.2.4), the 


most general solution for a(x,t) and e(x,t) for fixed Q must be of the form 


a(x,t) = ue + ver | (3.2.12) 


iNt-iqx 


e(x,t)=u,e +Vv emt (3.2.13) 


Substituting solutions (3.2.12) - (3.2.13) into (3.2.3) - (3.2.4) and using the 


fact that the modes are linearly independent, we obtain 


(-2+@h(k)q)u, +0 (k)qu, = 0 (3.2.14) 


(-2+04(k)p)v, +@(k)pv, = 0. (3.2.15) 


Using (3.2.10) and (3.2.11) we obtain 








ey ee (3.2.14a) 
2 9 &o 

ae Pe (3.2.15a) 
2 Mo Eo 


Equations (3.2.14a) - (3.2.15a) determine the ratios of u, to u, and v, to v, 
respectively. The ratio of u, to v, and u, to v. are determined by boundary 


conditions on the signal. 
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3. AMPLITUDE MODULTION TO FREQUENCY MODULATION 


Consider an arbitrary amplitude modulated signal at x=0 of the form 
a(0,t) = A[1+mcos(0t)|cos(at) , (3.3.1) 


where A is the amplitude, m is the modulation amplitude, Q is the modulation 
frequency and w is the carrier frequency. At the boundary, Eq. (3.2.13) 


becomes 
e(0,t)=u,e" + ve . 


Also, at x=0, the general expression for the phase becomes 


0(0,t) = te + Ste — at 
q p 


(x,t) 


08 
where we have integrated the expression, ae k, + a(x,t) for the modulated 


wavenumber and taken the associated arbitrary function of time in the result to 
be -wt. 
Because at x=0 the amplitude modulated signal is of the form (3.3.1), it 


follows from Eq. (3.3.3) 
inte 4 jet <9. 


This requires 





U,V 


—+4-1=0. (3.3.5) 
q Pp 
Because the nonlinearity is weak, p is approximately equal to q and 


u,=-V,. (3.3.6) 


which expresses the fact that energy is equipartitioned between the two 
propagating modes. 
At any point x, the energy «(x,t) is 


e(x,t) = U2 [cos(t - qx) + cos(Mt - px)] , (3.3.7) 


where we have considered the real part only. Using trigonometric identities and 


defining A = = and n= a , we obtain 
€(x,t) = 2u, cos Ax cos(Qt~ nx) . (3.3.8) 


If the modulation amplitude m is small compared to unity, from Eq. (3.3.7) 
£(0,t) = 2A2mcosnt . (3.3.9) 


Upon comparison with Eq. (2.5.2) we obtain 


u, =A’m. (3.3.10) 


Similarly, at any point x, the phase 6(x, t) is 
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0(x.t) = koX~ at - “= sin( Ax) e0s(2t~ nx) : (3.3.11) 





and from Eq. (3.2.14a) u, = 02 A2m. 
ole 


Also, using the fact that q = Q/«,’, we obtain the following expression for the 


amplitude modulated wave: 


a(x,t) = A[1+ 2mcos(Ax)cos(Qt- yx)] 


2 2 
X cos ot 2A esl foe sin(Ax)cos(Qt— nx) | . (3.3.12) 
O5e, 


Q 


Thus a consequence of the nonlinearity on an amplitude modulated signal is to 
introduce frequency modulation. 

Equation (3.3.12) exhibits reciprocity regarding frequency modulation to 
amplitude modulation. When cos Ax = 0 the signal is entirely frequency 
modulated. If we consider this point as the source of the modulated wave, it must 
follow that a consequence of nonlinearity on an frequency modulated signal is to 


introduce an amplitude modulation into the signal. 
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CHAPTER 4. SOUND WAVES IN A CYLINDRICAL WAVEGUIDE 


In this chapter we review the theory of sound wave propagation in a 
cylindrical waveguide. Following Morse and Ingard (1968), we examine the 
general linear theory and derive the waveguide dispersion relation. We discuss 
the perturbation technique of Keller and Millman (1972) that can be used to 
determine the effect of small finite amplitude of an acoustic wave. We specialize 
the results of the perturbation method to determine the nonlinear coefficient w » 
for the specific case of the first non-axially symmetric mode in a cylindrical 


waveguide. Lastly, we examine the experimental consequences of our results. 


1. LINEAR THEORY 


it 


The wave equation for simple harmonic waves p(x,y,z)e" propagating in an 


acoustic duct is 
a @ 8 o\ 
[+53 p-(2) p=0, (4.1.1) 


where p is the pressure, c is the thermodynamic speed of sound in the medium, 
and @ is the angular frequency of the wave. Setting the z axis parallel to the axis 


of the duct we seek a solution of the form 
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p='P(xy)e" , (4.1.2) 


where the factor 'Y satisfies the equation 


ax? * ay? Yr+y vY=0, (4.1.3) 
with 

ms 2 

(2) =k? +y?, (4.1.4) 


the waveguide dispersion relation for « (k). 
For a rigid wall duct, the fluid velocity normal to the surface wall must be zero. 


From the linearized momentum equation 
éu 
——=-Vp, 4.1.5 
a p (4.1.5) 


a harmonic wave satisfies 
ipo = Vp (4.1.6) 


where p is the density and u is the fluid velocity. Thus, in terms of the acoustic 
pressure the appropriate boundary condition for a rigid wall is 


“=0, (4.1.7) 


for every point on the wall perimeter with n the normal to the wall. 





There exists a discrete set of eigenvalues y, for which Eqs. (4.1.3) and (4.1.7) 
are Satisfied. From Eq. (4.1.4) it can be seen that to each mode corresponds a 
specific cutoff frequency cy ee For frequencies w greater than cy, k is a real 
number and the wave propagates along the tube. For frequencies less than cy,, 
k is imaginary and the mode does not propagate. 

The lowest eigenvalue 7, corresponds to the fundamental plane wave mode. 
For this mode, the phase velocity is the speed of sound in the medium. At 
frequencies below the cutoff frequency of the first higher mode only the 
fundamental plane wave mode propagates without attenuation. 

As the frequency is increased above the cutoff frequency of the first higher 
mode, the plane wave mode continues to propagate together with the first higher 
mode. 


The phase velocity of the jth mode is 


Ye eee (4.1.9) 


2 
Vo=—p=c i-{%) == (4.1.10) 


This velocity is less than c for frequencies above cutoff. 
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For a waveguide of circular cross section of radius R, the transverse 


eigenfunctions and eigenvalues of Eq. (4.1.3) are (Morse and Ingard, 1968) 


cos mwr 
= — |, 4.1.11 
raat (mb)dq 2 (4.1.11) 
TW 
Ymn = R’ (4.1.12) 


where r and 9 are the polar coordinates about the central axis of the cylinder 
and where each pair of subscripts m, n defines a wave mode propagating in the 
waveguide. The subscript m indicates the number of nodal planes while the 
subscript n indicates the number of nodal cylinders. The alternative possibilities 
of using cosine or sine functions express the fact that the non axially symmetric 
modes are degenerate. 

For each integer value of m there exist a set of values w which satisfy the 


boundary condition 
J! (nw) = 0. (4.1.13) 


The solutions w which satisfy this condition are denoted as w= a,,,. The plane 
wave mode corresponds to a9) =0, Jo (ma99)=1. From Eq. (4.1.12), the cutoff 
frequency of the (m,n)th mode is o,,,¢/2R Hz. 

The first mode above the plane wave mode is the (1,0) mode. This mode has 
one nodal plane extending radially out from the axis and no nodal circles. Since 
the axial velocity distribution is proportional to J,(xa, )/ R)cos(), the magnitude 


of the axial velocity along a line transverse to the nodal plane is 0 at the node 
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( d=+n/2) and a maximum at r=R. However, the axial velocities in the two halves 
of the waveguide are opposite in phase. This dipole mode is the dominant mode 
excited when the cylindrical waveguide is driven above cutoff with two speakers 
in an anti-symmetric push-pull configuration. 

Generally the wavenumber k may have real and imaginary parts and we can 
write 


4 
k === +inn, , (4.1.14) 
rj 


where i is the wavelength along the axis of the waveguide of the jth mode and 
n; is an attenuation factor. The sources of the attenuation are viscous and 
thermal losses at the duct walls as well as losses due to the acoustic admittance 
of the boundary surface. For a plane wave mode the attenuation due to 


thermoviscous losses is (Landau and Lifshits, 1959) 


= ve (14 JPr(e-1)) (4.1.15) 
c 


where v is the kinematic viscosity coefficient, Pr is the Prantl number, and & is 
the ratio of specific heats. An expression for the attenuation of the higher order 


propagating mode (1,0) as measured along the z axis is 


VOV Cp 
Ss SiO Prey; 4.1.16 
Mh = oe dli- ima r(&-1) (4.1.16) 
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The modification of the term above from its form in (4.1.15) can be understood 
by recognizing that the distance traveled by the wavefronts in the direction of 
their wavevectors is (c,/c) times greater than the distance measured along the z 
axis. Because of the circular geometry the wavefronts are not straight. This 


effect is : presented by the factor [1-(1/na.,,)*] . Thus, the viscous losses are 


multiplied by (c,/e[1-(1/ma149)?)). The contribution due to thermal losses remains 


the same because they are independent of the shear component of the velocity. 
The distance at which the pressure field decreases by a factor e' thee 

folding distance, is the inverse of the expression (4.1.16). It approaches 0 as the 

frequency approaches the mode's cutoff frequency, it has a maximum when 
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@ / Cy; 9 is approximately equal to 3'”, and goes as w’'? for high frequencies. 


2. THE NONLINEAR COEFFICIENT 

The combined effects of nonlinearity and dispersion may force a single 
wavetrain to split into two separate wavetrains moving with unequal velocities. A 
calculation of the magnitude of the nonlinear coefficient is thus required to 
determine whether the effect is measurable for the specific case of the (1,0) 
cylindrical waveguide mode. 

A general expression for the nonlinear coefficient for the wavenumber k, has 
been obtained by Keller and Millman (1970). They consider the effect of a small 
finite amplitude sound wave propagating in a rigid cylindrical waveguide. A 


perturbation method is used by expanding the velocity potential © and the 








pressure p in powers of the dimensionless amplitude « to solve the equations of 
fluid dynamics. 


Specifically, consider the Bernoulli equation 


aD 4 2 | dp 
—+=(Vo)°+ | —~=0, (4.2.1) 
at 2 ! p(p) 


and the continuity equation 
2 a 
p,[p(p)]V 7 — = avo. 2). 3(vo.v)(vo) . (4.2.2) 


where © is the velocity potential, p is the pressure p(p) is the density, p'° is a 


constant and p, is the derivative of the pressure with respect to the density. 


Equations (4.2.1) and (4.2.2) must be supplemented with the boundary 


condition 
—=0, (4.2.3) 


corresponding to the vanishing of the normal component of the velocity at the 
waveguide perimeter. 
Regarding the frequency w as given and k as unknown and taking the z axis 


as parallel to the waveguide axis, solutions of the form 


® = (x, y,kz- ot), p=p(x,y,kz— ot), (4.2.4) 


are sought. 
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Introducing the variable change z’ = kz and t’ = wt and then dropping the 
g 


primes from the resulting equations allows us to rewrite Eqs. (4.2.1) and (4.2.2) 


in the form 


eD | ao 0° 
—— +k* —~ ]-o 
polo] So+ So 42 S$ 
V5 9, 92 & OG 2M @\ ro Ho |, do 
(2° at ax ax ey by éz 6z }\ dx? dy? éz2 ) 
where 
o(x,y,z-t+2x) = 0(x,y,z-t), 


p(x, y,z—t+2n)=p(x,y,z-t). (4.2.7) 


The Taylor expansion of ®, p, and k in terms of ¢ can be obtained as follows 


(4.2.8) 


() 8 
®,(%y,Z-te) =e, +FOn te 


P,a(% y,Z-te) =p? +ep"” + (4.2.9) 


(1) 


k,(e) =k) + ek! 2)i (4.2.10) 


2 
€ 
vat ested -§ 


where the superscripts denote the jth derivative with respect to « and the 
subscript denotes the propagating mode. The justification for Eq. (4.2.10) is that 


the Taylor expansions of ® and p in terms of « do not exist unless k also 


depends one. 








Keller and Millman (1970) insert Eqs. (4.2.8) - (4.2.10) into Eqs. (4.2.5) and 


(4.2.6). They show that k”) = 0 and obtain the following general expression for 


the coefficient K?) 


(2) 1 
kK, =~ [(oV 7B Vit, 
(0) n n 

2c*k, 


3 0) ]* 
~ 3 YY Yn 4 oh} ea IVF 
2 
20K. + 2m} +P) dxdy . (4.2.11) 


(0) 


where c is the speed of sound in the medium, k, 


is the linear wavenumber for 


2 2 
which [k |" = oy? ( 7, is the eigenvalue of the wave equation (4.1.3) and 
n C2 n 


Y,, is the corresponding eigenfunction), V; is the transverse Laplacian, B, is 


determined by solutions to 


(v2 +4? )p, 2 (2) {Key +2n)¥2| (4.2.12) 


and 














2 2 
| aw, ota Mn ay, (2) ary, (4.2.13) 


Wa Ya Yaa = [ OX ax? ax ay axdy ay ay? 


For an ideal gas equation of state p(p)=Ap* 
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2 
h=(—-1)— =, 4.2.14 
6-175 (4.2.14) 


eg _ = 4 eS 2 2 
§2 140) eee mes Dar g Seat e {(ee P28}. (4.2.15) 


We now specialize these results to the case of the first non-axially symmetric 
mode propagating in a cylindrical waveguide of radius R. From Eq. (4.1.11) the 


normalized solution of the full wave equation for this mode is 





R 


RJ, (org )af(moryg)?-1 


This mode is doubly-degenerate and even when an external drive breaks the 


V270 49 COS()u, | 


symmetry, nonlinearities resonantly couple the two eigenmodes. However, 
because a real waveguide is not perfectly circular the degeneracy is removed 
and the frequency difference is a monotonic function of the eccentricity. Thus 
we will assume a nondegenerate solution of the form (4.2.16) and determine the 
nonlinear coefficient (4.2.11) accordingly. 

To obtain an expression for B, . we insert (4.2.16) into (4.2.12) and transform 
to polar coordinates. The right hand side of the inhomogeneous Eq. (4.2.12) is 


of the form F(r) + H(r) cos(2o), where F(r) and H(r) are in terms of products of 
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two Bessel functions. The solution of (4.2.12) satisfying the boundary condition 
(4.2.3), which becomes B, o(R) = 0, is 


R R 
B,o(r) = [Go (r,s)F(s)sds + cos(2$){ G, (r,s)H(s)sds : (4.2.17) 
0 0 


Here, the Green's functions G,,(r,s) are given by 


mance Jin (270. of. /R) 
me 2 dy" (2 70010) (4.2.18) 


X [Jin (270 of PRYNmy (270.10) —Non (27004 of, HR) (220410 )] . 


where N,, is the Neumann function of order m, and r, (r.) stands for the greater 


(smaller) of rand s. Using (4.2.14)-(4.2.18) we obtain for k\? the following 


4 
k?(@) = V2n0 i 
ve RJ, (09 Jaf (HL, 9)” -] 


2 2 

T A19 
2p 21-(0) 
2c°R Kio 


(4.2.19) 


{0.305 +.0.38x? - 1.463x* +0.371x5} 


where x = @R/(cza, 9) is the ratio of the frequency to the cutoff frequency. The 
expression (4.2.19) has meaning for x>1. It is negative for 1< x < 1.897 and 
becomes positive for frequencies larger than 1.897 times bigger than the cutoff 


frequency. 





The wavenumber is the result of developments in powers of « which is the 
amplitude of the linearized solution. In terms of the ratio of the acoustic pressure 


to the ambient pressure we obtain 





” {z ah)" atrial (a1) (O49)? - ' a(2 ). (4.2.20) 
c R 2 V270 49 Eo Po 
The propagation constant k, , depends upon amplitude. It is less than the linear 
value for frequencies between the cutoff frequency and 1.897 times the cutoff 
frequency. According to the physical arguments presented in the Chapter 1, a 
modulation of a wavetrain will be stable and an intial disturbance will split into 
two disturbances moving with different group velocities. On the other hand, for 
frequencies greater that 1.897 times the cutoff frequency the phase velocity is a 
decreasing function of amplitude and modulations of a wavetrain will become 
unstable to the formation of solitons. 

The formalism developed in the previous chapters relies on the determination 
of the nonlinear coefficient w,. Equation (4.2.20) gives the fequency 
dependence of the propagation constant k,,. Thus, we can solve implicitly for 
the frequency as a function of wavenumber k = k,,. If we set @=Oyt2(P/P,)*, it 


can be shown that 


onto =o S| eeehleear 


RJ, Saieeeideeal >) k2(o se 2) , (4.2.21) 
0 


V2 M049 &@o(k) 
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where 





ra : 
@ (k)=¢ “| +k? (4.2.22) 


is the linear dispersion relation and Wo! is the linear group velocity at the 


wavenumber k. 


3. EXPERIMENTAL CONSEQUENCES 

For a fixed frequency, Eq. (4.2.20) predicts that the wavelength of a nonlinear 
wave changes with amplitude. It increases with amplitude for frequencies less 
than 1.897 times the cutoff frequency and decreases with amplitude for 
frequencies that are larger than 1.897 times the cutoff frequency. The mode is 
modulationally stable in the former case and we will restrict our analysis to this 
frequency range. 

The effect of increasing the wavelength with amplitude can be determined 
experimentally by measuring the relative phase between two points. ff the 


distance d between two points is less that the wavelength 


= (4.3.1) 
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where Aq is the change in phase difference @ between two points due to a 


change in amplitude, and Ak =k,, —k‘0. The effect can be magnified if the 


distance between the two points several wavelenghts apart. In this case 
Ao Ak 


a 4.3.2 
o+2nN_ k ( 





where N is the integer number of wavelengths in the distance d between the two 
points. The slope of the measured relative phase between two points as 
function c*t the square of the amplitude will thus yield an experimental value of 
the nonlinear coefficient. This value in turn, can be used to determine the 
separation in time as a function of distance of the split pulses. 

To determine the difference of the times of arrival at a given location of the 
two disturbances that result from the splitting of an initial disturbance we can use 
Eq. (3.1.4). In terms of the ratio p/p, of the acoustic pressure to the ambiert 


pressure, (3.1.4) becomes 


W/2 
G2 _ oy (K)#{0, (kag! (K)} | B (4.3.3) 
dt Po 
A finite extent modulation of a nonlinear wavetrain generated at a source at z=0 
will split into two disturbances that arrive at a distance z with a time difference 


given by 


i“ 1/2 
At, (2) = z2ioe oa) ( P (4.3.4) 
0 
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In order to determine if this effect is observable, we must also determine if the 
time spread due to dispersion is smaller than the time due to splitting. 

The average velocity of a smooth pulse or a finite wavetrain that does not cut 
off too violently is the group velocity Vg=dwo/dk=ao’ . Due to dispersion the pulse 
spreads because different frequency components move with different phase 
velocities. The spreading of the pulse can be accounted for by noting that a 
pulse with initial temporal width At, must have an inherent spread of frequencies 


Aw ~ (1/Aty). This means that the group velocity spread is of the order 


/ 
Av, ~ Aw ~ 22 (4.3.5) 





which at a point z implies a spread in time of the order of Av, z / (wo')*. Because 
the uncertainty in time is the square root of the sum of squares the width At(z) at 


location z is 


(4.3.6) 





The expression for At(z) shows that in a medium with dispersion a narrow pulse 
spreads rapidly because of its broad spectrum in frequencies. Thus the criterion 
for a small change in shape is that @/, At, >> o/!. 


Thus, for the splitting to be observable for al! distances we require that 
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2atyo Pes), Ears . (4.3.7) 
@ Po 
For a 10 ms disturbance superimposed on a signal with a frequency 10% higher 
than the cutoff frequency the minimum ratio of acoustic pressure to ambient 
pressure is 0.005, or 148 dB (re 20 pPa). 

The estimate (4.3.7) guaranties that the curves (4.3.4) and (4.3.6) do not 
intersect. For low intensity levels the two curves intersect. However, it should 
still be possible to observe the splitting of a pulse into two pulses for distances 


less than the distance given by the real solution to 


1 Hee 4p 3 
tf cs | = [at 22les al" | , (4.3.8) 


23 


9 Aty (@)? (Po 


provided the value At, given by (4.3.4) is within the limits of instrumental 
resolution. Spreading due to dispersion will dominate for distances larger than 
the distance determined by the real solution of (4.3.6). However, even when the 
splitting of a disturbance into two cannot be observed, the effects of the double 
group velocity will be to spread a given pulse beyond the values predicted by 


linear dispersion. 





CHAPTER 5. EXPERIMENTAL RESULTS 


In this chapter we discuss the experimental results related to propagation of 
the (1,0) cylindrical waveguide mode in the frequency regime for which the mode 
is modulationally stable. We briefly describe the experimental apparatus used in 
our experimental investigations. The experimental results include 
measurements of the (1,0) mode linear cut off frequency, group velocity, and 
attenuation constants for frequencies above the cutoff frequency. From 
frequency response measurements we also determine the relative amplitudes of 
the plane wave and dipole modes. Also, we measure the spreading of a low 
amplitude acoustic pulse. Finally, investigations of nonlinear effects including 


the determination of the nonlinear coefficient at 4200 Hz and preliminary time 


domain observations of group velocity splitting in amplitude modulated and in 


pulse modulated acoustic waves are presented. 


1. EXPERIMENTAL APPARATUS 

The experimental apparatus is a traveling wave tube made of an extruded 
aluminum pipe of average inner diameter 1.957 inches, outer diameter 2.377 
inches, length 70 feet, which consist of seven 10 foot sections joined by collars 
and flanges. An anechoic end termination made of steel wool with a tapered 


density distribution spans the length of the 7th section. Absorber and junction 





reflection amplitudes are typically -40dB down for frequencies between 1kHz and 
4kHz (Dorff, 1991). 

Each of the first six sections of pipe have 4 microphone ports which 
accommodate ENDEVCO series 8510B piezoresistive pressure transducers. 

The ports have a uniform spacing of 2.5 feet. An additional port was added 
directly below the 2nd port in the 5th section (from the drivers’ end of the pipe). 
This port adds the capability of instrumentally eliminating the detection of the 
plane wave mode by connecting the microphones with opposite polarity and 
adding their output. 

Two JBL model 2445H compression drivers were mounted at one end of the 
wave tube via a two-to-one adapter. The first microphone port is at a distance of 
0.83 m from the drivers. The drivers are titanium diaphragm drivers rated at 100 
watts from 500 Hz to 1000 Hz and 150 watts above 1 kHz. 

A two-to-one adapter was designed such that the drivers could be mounted 
directly opposed (Fig. 5.1). The inner diameter of the adapter was machined to a 
diameter of 1.975 inches to match the inner end diameter of the wave tube. The 
speakers were mounted flat to the adapter onto rectangular attachment flanges. 
A 2 inch to 1 inch tapered bore made through the attachment flanges and the 
adapter connected the 2 inch driver opening to the 1.975 inch diameter 
waveguide opening. 

In this configuration, the drivers operate in push - pull when driven 180 
degrees out of phase. The plane wave mode can in principle be nulled by 
adjusting the amplified signal to each driver. In this manner, at frequencies 


above cutoff the (1,0) mode is the main propagating mode. 
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Fig 5.1 Traveling Wave Tube 
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2. EXPERIMENTAL RESULTS 
a) Frequency Response 

In an acoustic waveguide the plane wave mode propagates at all 
frequencies. Although the experimental procedure calls for nulling the plane 
wave component of the acoustic field, the frequency responses of the drivers are 
not identical. Therefore, even if the plane wave component is nulled at one 
frequency, it may not be so for other frequencies in the sweep. 

The role played by the plane wave field was investigated by taking 
frequency response measurements with a swept sine source at the first port 0.83 
m from the drivers’ end and at each alternate port. The frequency span of the 
sweep is from 3900 to 4400 Hz, with a resolution of 401 points per sweep. The 
plane wave component was nulled at 3900 Hz. 

To investigate the effects that the presence of a plane wave might have in 
the dipole field, we obtained frequency responses measured at the double ports 
located at 13.79 m from the drivers' end. The double port arrangement allows 
measurements to be obtained by the output of a single microphone as well as 
the output of the sum of the two microphones with the same and with opposite 
polarities. 

The frequency response measured from a single microphone is the 
superposition of a plane wave and a dipole field. For a given frequency, the 
corresponding wavelengths of the two modes are not the same. Consequently, 
the relative phase of the corresponding acoustic pressure is a function of 


frequency. Thus we should expect that the interference between the two fields 








at a given point be reflected in the frequency response at this point. Fig. 5.2a 
shows a sinusoidal variation of the amplitude as a function of frequency. 

Fig. 5.2b shows the output at the dual port of two microphones with the 
same polarity. In this arrangement the dipole field is canceled and the output 
corresponds to a frequency response of the plane wave. The inconsistency in 
the amplitudes for the frequency responses of Fig. 5.a as compared to that in 
Figs. 5.b,c is due to the difference in the microphone gain used to obtain these 
signals. The microphone gain in Fig. 5.a is 500 while the gain used to obtain 
Figs. 5.b,c is 10k. Note that for frequencies above 4050 Hz, the amplitude 
decreases with increasing frequency. In Figs. 5.2a,b the frequency response 
abruptly dips at 4030.2 Hz which corresponds to the cutoff frequency for the 
(1,0) mode. 

Alternatively, Fig. 5.2c shows the output at the dual port of two 
microphones with opposite polarity. In this arrangement the plane wave field is 
canceled and the output corresponds to the frequency response of the (1,0) 
mode dipole field. For frequencies greater than 4200 Hz, the frequency 
response is essentially flat. 

A comparison of Figs. 5.2a,b&c leads us to the conclusion that 
interference is occurring between the plane wave and the (1,0) mode. As shown 
in Figs. 5.2b&c the plane wave mode is observed to be between 9 and 25 dB 


down from the dipole mode for frequencies above the cutoff. 


47 


{dBVrms) 


Amplitude 


Amplitude (d8Vrms) 


Amplitude (dBVrms) 


Fig. 5.2 


f 
ae OTT 
VE VS AV AR ia a 


Wika 
Ne 





[i 4200 
Frequency (H2z} 


Frequency response of the traveling wave tube measured at 
the dual port 13.79 m from the drivers. a) Output from 
single microphone. Mike gain is 500. b) Output from two 
microphones with equal polarity c) Output from two 
microphones with opposite polarity. Mike gain for b,c is 10k. 
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b) Attenuation Constant 

From the frequency response at each port we obtained microphone signal 
amplitudes for 401 frequencies between 3900 Hz and 4400 Hz. A matrix was 
constructed from this data. The rows of the matrix contain the signa! amplitudes 
at fixed frequency for 14 microphone port locations. The columns list signal 
amplitudes at a fixed location for 401 separate frequencies between 3900 and 
4400 Hz. Using this matrix, we obtained plots of the amplitude versus distance 
for each frequency. From an exponential fit performed on the data, we obtained 
the attenuation constants for frequencies between 4050 and 4400 Hz. 

Fig. 5.3a is a plot of the e"' folding distance as a function of frequency. 
The attenuation constant corresponds to the inverse of this distance. The plot of 
this data shows irregularly high and low e"' folding distances between the 
frequencies 4100 to 4200 Hz. As suggested by Fig. 5.2 the variation in the e” 
folding distance between these frequencies is due to the rapid variations in the 
frequency response at the different ports between these frequencies. That is, 
the frequency corresponding to a peak at one port may not necessarily 
correspond to a peak at another port. This observation is supported by the fact 
that variation in the e"' folding distances is not observed at frequencies where 
the amplitude of the frequency response is relatively flat. 

In Fig. 5.3b we have removed the extreme high and low attenuation 
factors in the frequency range 4100 to 4200 Hz. The solid curve in Fig. 5.3b, 
obtained from Eq. (4.1.15), corresponds to the predicted attenuation length due 
to thermoviscous losses . Excellent agreement between theory and experiment 
is observed up to 4300 Hz where pipe radiation losses possibly become 


important. 


49 


1le folding distance (m ) 





400 4200 4300 
Frequency ( Hz ) 


1le folding distance (m ) 





4050 410 48 4200 4250 4300 4350 


Fig. 5.3 Plot of the 1/e folding distance as a function of frequency. 
The extreme high and low variations in (a) have been 
removed in (b). The solid line corresponds to theory. 
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c) Group Velocity Measurements 
The frequency response plots of Fig. 5.2 show a change in behavior for 
frequencies above the cutoff frequency of 4030.2 Hz. According to the theory, 
when the drivers are in a push-pull configuration the dominant propagating mode 
is the (1,0) mode. For this mode, the group velocity depends on frequency 


according to the relation (4.1.10) 
5 \1/2 
v,=¢(1-(o,/0)*) (5.2.1) 


where a, is the cut off frequency of the (1,0) mode. Thus for a fixed distance D, 


the time of travel At is 
(1/at)? =¢?/D? {4-0,2(Yo?)} (5.2.2) 


We made measurements of the time of flight of a 100 cycle burst with 
carrier frequencies above cutoff. The envelope was sinusoidal. The time delay 
measurements were made at a single port and a fixed location 13.79 m from the 
drivers on a Nicolet Oscilloscope triggered externally from the driver signal. 

Time delays were measured from the time of trigger to the time of half maximum 
signal amplitude for five frequencies (Coppens, 1966). Fig. 5.4 shows the plot of 
the resulting data. As predicted from (5.2.2) there is a linear relationship 
between the inverse of the square of the time delay ( 1/At?) and the inverse of 
the square of the frequency ( 1/f* ). A least square fit of the data yields a cutoff 


frequency of 4040 Hz. This result is within 0.2% of the value obtained from the 


frequency response data. The experimentally obtained value for the speed of 


sound is 338.7 m/s. 





1iAt? (sec? ) 





Fig. 5.4 Group velocity measurements. The graph shows a linear 
dependence between the inverse of the square of the time 
delay and the inverse of the square of the frequency for the 
(1,0) mode for 5 frequencies above cutoff. Slope = 877.9 ; 
Intercept = -1.43 x 107°. 
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d) Measurements of Pulse Spreading 

Equation (4.3.6) estimates the spreading of an acoustic pulse due to 
dispersion. This equation predicts that the spreading is significant only at 
distances much greater than the length of our waveguide or for frequencies very 
close to cutoff. Specifically, predictions for the spreading of a 21.7 ms pulse 
propagating a distance of 7.69 m range from 0.002 ms for a 4.6 kHz signal to 0.1 
ms for a 4.2 kHz signal. These times are within the instrumental uncertainty of 
the Nicolet oscilloscope. 

in order to verify the prediction that the spreading due to dispersion was 
negligible, we analized the speading of a low amplitude N cycle pulse of 21.7 ms 
duration. The carrier was modulated to obtain a smooth sinusoid envelope. The 
number of cycies N was varied to keep the input pulse duration constant. 
Measurements of pulse spreading were made by measuring the total time 
duration of the pulse signals obtained at locations 0.83m and 7.69m from the 


drivers. Table 5.1 shows the measurements obtained for 5 frequencies. 


Table 5.1 
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The data verify that for the propagation distances and frequencies 
measured the contribution to the spreading from dispersion is negligible. 

The spreading of the pulse at 0.83 meters is apparently due to the fact 
that at this distance the plane wave has not yet unraveled from the center of the 
dipole mode pulse. That is, the microphone signal consists of a plane wave 
component as well as a dipole component where each travel with different 
velocities of propagation and therefore arrive at the microphone at different 
times. For example, in considering the data for the 4.4 kHz signal, the plane 
wave component arrives at the microphone in 2.45 ms. The dipole component 
arrives in 6.1 ms. The combined signal would therefore appear to be 3.6 ms 
longer than the input signal. Experimentally, the signal at the .83 meter location 
for 4.4 kHz is observed to be 3 ms longer than the input. 

Signal distortion prevented a reliable measurement of the duration of the 


4.2 kHz signal from the 7.69m location. 
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e) Measurement of the Nonlinear Coefficient 

As shown in Fig. 5.4, for frequencies above cutoff the mode is dispersive. 
For small but not infinitesimal amplitude, the dispersion relation for the (1,0) 
mode is of the form given by Eq. (2.2.1). According to this expression, for fixed 
frequency a change in the amplitude induces a change in wavelength of the 
mode. If the nonlinear coefficient is positive the changes in the wavenumber 
decrease linearly with the square of the amplitude. Thus, by measuring the 
relative phase between two points as a function of amplitude we can determine 
both the validity of Eq. (2.2.1) and the nonlinear coefficient. 

Using a PAR 5210 lock-in amplifier we determined the amplitude of the 
acoustic field at 13.79 m from the drivers and the phase of this signal relative to 
a reference signal located at a distance of 9.98 m from the drivers. By varying 
the driver amplitude at a fixed frequency, in this case 4200 Hz, we obtained the 
relative phase between the two signals as a function of the amplitude squared. 
The plane wave component of the acoustic field was nulled at this frequency. 

The data in Fig. 5.5 are in terms of the acoustic pressure at the source. A 
reasonable estimate of the amplitude at the source was obtained by multiplying 
the calibrated acoustic amplitude at 13.79 m by the inverse attenuation factor 
e'9:79/339 where 3.39 m is the 1/e folding distance of a (1,0) mode of frequency 
equal to 4200 Hz. 

Fig. 5.5 shows a linear relation between the relative phase and the square 
of the acoustic amplitude at the source. In order to get an idea of the intensity 
levels used for this measurement, the pressure amplitude has been normalized 
to the atmospheric pressure. The plotted range varies between 150 dB and 170 


dB (re 20 »pPa). The nonlinear coefficient can be determined from using the 
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slope of the line and Eq. (4.3.2). Fora frequency of 4200 Hz, the theoretical 
value predicted by Eq. (4.2.20) yields Ak = -12.58 (p/p,)*. In the span of 3.81 m, 
corresponding to the distance between the two ports there are 13 wavelegths. 


Thus we obtain 


2 
BK +2nN) = -1.0} 2.| ; (5.2.3) 
k Po 


The predicted slope is 12% bigger than the experimental value of - 41.36. 
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Fig. 5.5 Plot of the relative phase versus the square of the 
normalized acoustic pressure at the source. Pressure 
amplitudes are normalized to 1 atmosphere. 
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f) Measurements of Group Velocity Splitting 

An initial disturbance superimposed on a signal with a frequency 10% 
higher than the cutoff frequency splits into two and the differential time of arrival 
is approximately 30 z(p/p,) ms at a distance z from the source (Eq. (4.3.4)). 
Fig. 5.6 shows the microphone output on a Nicolet digital oscilloscope for three 
different locations and three different intensity levels. 

The input signal is generated using 2 HP 3314 function generators. A 
4500 Hz carrier is amplitude modulated by a 100 Hz inverted half-sine wave at 1 
sec intervals. The first, second and third rows in Fig. 5.6 correspond to intensity 
levels of 140 dB, 144 dB, and 156 dB at the source, respectively. At 140 dB no 
apparent splitting of the negative pulse is observed up to the shown distance of 
7.69, while the spreading of the pulse is noticeable. However at 144 dB and 156 
dB the pulse undergoes spreading and a splitting becomes readily apparent for 
distances of 4.64 m and beyond. 

We performed quantitative analysis of the signals displayed in Fig. 5.6 
The procedure used consists of multiplying the data by the cosine of the carrier 
frequency, taking the Fast Fourier Transform (FFT) of the product and then 
filtering the modulus of the FFT to remove the high frequency components of the 
frequency spectrum. The inverse FFT of the filtered frequency spectrum is then 
taken to recover the signal modulation exclusively. 

With this procedure we can determine the time spread of the pulse and 
the differential time of the split portions. Measurements were made from the 
waveform peaks. For a sound level xf 140 dB a 10 ms pulse spreads into an 


11.8 ms pulse at 4.64 m and into a 14.2 ms pulse at 7.69 m. Thus the delayed 
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time is 2.4 ms and it is not at all described by dispersion. Furthermore, at 10.73 
m the pulse undergoes splitting where the difference between the two peaks is 
5.5 ms. 

Similarly, for an intensity of 144 dB the apparent spread of a 10 ms pulse 
at 4.64 m is 14.7 ms, while at 7.69 m the spread is 17.0 ms. This corresponds to 
a delay of 2.3 ms. The split pulses at these locations are 4.1 ms and 6.1 ms 
apart respectively. Remarkably there is a delay of 2 ms in the split pulses 
between these two locations. That is, the splitting the pulse occurs with a delay 
nearly equal to the differential spread. 

Furthermore, a comparison of the 144 dB and the 156 dB signals at 7.69 
m shows a significant increase in the time interval between two peaks of the 
modulation. Fig. 5.7 shows that for the lower amplitude signal, the time interval 
between the two most negative peaks of the modulation is 6.1 ms. For high 
amplitude the time interval between these same two peaks is 8.4 ms. The time 
delay is 2.3 ms verifying that the splitting increases as a function of amplitude. 
On the other hand, the 17 ms spread is the same for both amplitudes. 

Using a single HP 3314 function generator in the N cycle mode, we 
generated a pulsed sine wave of 50 cycles with frequency 4.2 kHz. We varied 
the amplitude of the driver input signal and obtained time domain measurements 
of a microphone output signal at 10.74 m from the drivers. The envelope of the 
pulse was recovered using a procedure identical to that described above. 

Fig. 5.8 shows the modulations of a 150 dB and a 178 dB bursts. The 
time intervals between the two marked peaks are 33.75 ms and 35.56 m 


respectively. The time delay is 1.8 ms. 
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Fig. 5.7 Modulation envelope of a 4500 Hz carrier wave at 7.69 m 
the source showing peak to peak time intervals. The 
intensity level at the source is: a) 144 dB, and b) 156 dB. 
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Fig. 5.8 Envelope a 50 cycle burst at a frequency 4.2 kHz 
at 10.74 m from the drivers. The acoustic intensity at 
the source is a) 150 dB and b) 178 dB. 
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CHAPTER 6. CONCLUSIONS, POTENTIAL APPLICATIONS, AND 
FUTURE WORK 


1. CONCLUSIONS 

This work reports the results of theoretical and experimental investigations 
involving the phenomena of group velocity splitting. We have presented a 
quantitative derivation of group velocity splitting and given a comprehensive 
descriptior. of the process whereby a modulation of finite extent splits into two 
separate disturbances traveling with two different velocities of propagation. Also, 
we have predicted that for modulationally stable nonlinear dispersive systems, a 
consequence of group velocity splitting is to introduce a frequency modulation 
component into an amplitude modulated signal. 

We have made the first experimental observations of the change in 
wavelength with amplitude for a mode above cutoff in an cylindrical acoustic 
duct. The experimental value obtained is within 12% of the theoretical 
prediction. 

To our knowledge, we have made the first experimental observation of pulse 
splitting in a nonlinear dispersive medium. Experiments show that a single 10 
ms duration high amplitude (156 dB ) pulse propagating 7.69 m down the 
waveguide splits into two pulses traveling at different speeds. Specifically, a 
single peaked modulated acoustic signal is observed to split into two peaks 


separated in time by 4.1 ms in fair agreement with theory. Additionally, since the 
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spreading of the pulse at these distances is not described by dispersion, we 
conclude that the observed spreading of the modulation at small propagation 


distances is also a result of group velocity splitting. 


2. POTENTIAL APPLICATIONS 

The importance of establishing experimentally the effect of group velocity 
splitting lies in the possible application to fiber optic communication, command, 
and contro} links. For example, potential applications could include use of pulse 
splitting in fiber optic networks requiring increased data redundancy. 

The linear group velocity of guided waves in dielectric fibers is determined by 
a combinatior. of the material dispersion and geometrical waveguide dispersion. 
In the frequency range of visible light and for a small difference between the clad 
and the core material, the material dispersion dominates and the dispersion is 
negative. However if one uses cladding having much smaller index of refraction 
than the core or longer wavelengths a positive dispersion results because 
geometrical dispersion dominates in the first case and material dispersion 
becomes anomalous in the second. Under these conditions pulse splitting 
becomes possible. 

For optical fibers with a positive dispersion, another potential application of 
the effects presented in this work is the possibility of an all optical AM-FM 


conversion. 
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Finally, the effects of group velocity splitting in coastal shallow water 
acoustics should be studied. For example, potential applications exist in the use 


of pulse splitting as a littoral environment mine countermeasure. 


3. FUTURE WORK 

Future work should include continued investigation of the pulse splitting 
phenomena using input signals of smooth modulation. The splitting appears to 
be more easily observable for pulses of this type. 

Future experiments should also include determination of the relationship 
between the relative phase and the square of the wave amplitude for a range of 
frequencies above cutoff. As predicted by the theory presented in Chapter 4, the 
nonlinear coefficient changes sign for frequencies above 1.897 times the cutoff 
frequency. For this frequency range, modulational instabilities lead to the 
formation of solitons. To our knowledge, no experimental observation of 
acoustic solitons has been made. A future experiment should explore this 
possibility. 

Experimental verification of the predicted AM-FM conversion effect should 


also be pursued. 
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